An extension of the finite difference time domain is applied to solve the Schrödinger equation. A systematic analysis of stability and convergence of this technique is carried out in this article. The numerical scheme used to solve the Schrödinger equation differs from the scheme found in electromagnetics. Also, the unit cell employed to model quantum devices is different from the Yee cell used by the electrical engineering community. A bound for the time step is derived to ensure stability. Several numerical experiments in quantum structures demonstrate the accuracy of a second order, comparable to the analysis of electromagnetic devices with the Yee cell.
展开▼
